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Abstract 

The isosinglet unpolarized and isovector polarized twist-2 quark distributions of the 
nucleon at low normalization point are calculated in the large— N c limit. The nucleon 
is described as a soliton of the effective chiral theory. We derive the expressions for the 
distribution functions in the large-iV c limit starting from their definition as numbers of 
partons carrying momentum fraction x in the infinite momentum frame. We develop 
a numerical method for computation of the quark and antiquark distributions as sums 
over the quark single-particle levels in the pion field of the soliton. The contributions 
of the discrete bound-state level as well as the Dirac continuum are taken into account. 
The quark- and antiquark distributions obtained explicitly satisfy all general require- 
ments. Results are in reasonable agreement with parametrizations of the data at low 
normalization point. 
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1 Introduction 



The evolution of parton distributions with q 2 in the asymptotic region is well understood to- 
day, being governed by the renormalization group equation of perturbative QCD. A complete 
description of experiments at large q 2 requires, however, the knowledge of parton distribu- 
tions in the nucleon at some initial normalization point. Several sets of input distributions 
were determined by fits to the experimental data at large q 2 [[TJ |2], ^ f2j . All these fits include 
antiquarks and gluons at a low normalization point. 

Recently, we have formulated an approach to calculate the twist-2 parton distributions at 
low normalization point in the limit of a large number of colors (N c ), where the nucleon is 
described as a chiral soliton ||. At low energies, QCD may be approximated by an effective 
theory whose degrees of freedom are quarks with a dynamically generated mass, interacting 
with pions, which appear as Goldstone bosons of the dynamically broken chiral symmetry. 
The nucleon emerges as a classical soliton of the pion field p. This picture is known to give a 
successful description of hadronic observables such as the nucleon mass, magnetic moments, 
form factors etc. [[/J. In we have shown that this approach possesses all necessary req- 
uisites for a successful description of the leading-twist parton distributions of the nucleon. 
The normalization point of the distribution functions obtained in this way is of the order 
of the ultraviolet cutoff of the effective chiral theory, typically ~ 600 MeV. Let us briefly 
summarize the main characteristics of this description ||: 

i) Classification of quark distributions in the large N c -limit. In the large- A c limit the 
quark distributions are concentrated at values of x ~ 1/N C . Combining this fact with the 
known large A^-behavior of the integrals of the distributions over x, one infers that the quark 
distributions in the large N c - limit can be divided in "large" and "small" ones. The leading 
distributions are the isosinglet unpolarized and isosinglet polarized distributions, which are 
of the form 

D l ^ c (x) ~ N 2 p(N c x), (1.1) 



where p(y) is a stable function in the large N c - 
bution considered. The isovector unpolarized 
only in the next-to-leading order of the 1/N C - 

D small (x) ~ 



limit, which depends on the particular distri- 
and isosinglet polarized distributions appear 
expansion, and are of the form 

N cP {N c x). (1.2) 



ii) Sum rules and antiquark distributions. The chiral soliton model is a field-theoretic 
description of the nucleon, which preserves all general requirements on parton distributions. 
In particular, the standard sum rules for parton distributions and their positivity properties 
are satisfied automatically within the model. Also, a consistent description of the antiquark 
distributions can be achieved in this approach. 

Hi) Parametric smallness of the gluon distribution. When working with the effective 
chiral theory, it is implied that the ratio of the dynamical quark mass, M, to the UV 
cutoff, A (not to be confused with the QCD scale parameter, Aqcd) ; is parametrically small. 
For M/A <C 1, the quark distributions computed in the effective theory may be identified 
with the "current" quark distributions of QCD. The gluon distribution is zero at this level, 
more precisely, it is 0(M 2 /A 2 ). For finite M/A, the quark distributions computed in the 
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effective theory should be interpreted as distributions of "constituent" quarks — objects 
which themselves have a substructure in terms of QCD partons. The gluon distribution 
inside these objects could in principle be recovered from the effective theory if one knew the 
precise way how the UV cutoff arises as a result of integration over the original QCD degrees 
of freedom. These statements can be made more precise in the framework of the instanton 
vacuum, which on one hand allows to derive the effective chiral theory, on the other hand 
can be used to evaluate the gluon distribution directly, using the method developed in [|J. 
One finds that the gluon distribution is suppressed relative to the quark distributions by a 
factor of the packing fraction of the instanton medium ||. 

In this paper, we study the properties of the iV c -leading quark and antiquark distribu- 
tions, namely the isosinglet unpolarized and isovector polarized distributions, in the approach 
formulated in ||. First, we rederive the basic formulas for the parton distributions in the 
effective chiral theory in a new way. In || these formulae were obtained from the exact 
QCD expressions for the parton distributions as matrix elements of quark bilinears with a 
light-like separation || [TIJ. In this paper we take the original Feynman point of view [IT 



that parton distributions are given by the number of partons carrying a fraction x of the nu- 
cleon momentum in the nucleon infinite-momentum frame. Despite the apparent difference 
in wording we show here that the two definitions are, in fact, equivalent and lead to identical 
working formulae for computing parton distributions. We think it is remarkable that the 
actual equivalence of the two well-known definitions can explicitly be demonstrated within 
this field-theoretical model of the nucleon. The deep reason for the equivalence is that the 
main hypothesis of the Feynman parton model, namely that partons transverse momenta do 
not grow with q 2 [pj| , is satisfied in the model under consideration. 

Second, we investigate the influence of the ultraviolet cutoff of the effective chiral theory 
on the distribution functions. This not only includes the asymptotic dependence on the 
cutoff parameter (i.e., the UV divergences), but also, and more importantly, the dependence 
on the regularization scheme adopted to make the distributions finite. It is crucial that the 
regularization method not violate the completeness of quark single-particle states in the 
background pion field, in order to preserve causality, that is, the anticommutation relations 
of quark fields at space-like separations. A regularization which meets this requirement is, 
for example, the Pauli-Villars subtraction. We show explicitly that it leads to quark and 
antiquark distribution functions satisfying all general requirements, such as rapid decrease 
for large x, uniform logarithmic dependence on the cutoff, and positivity. On the other 
hand, regularization methods based on an energy cutoff, such as the popular proper-time 
regularization of the determinant, violate causality and lead to unacceptable results for the 
distribution functions. 

Third, we develop a numerical method for exact computation of the quark and antiquark 
distribution functions as sums over quark single-particle levels in the background pion field. 
(In H the polarized distributions were estimated using an approximation, the so-called 
interpolation formula). Since the quark distribution functions are given by matrix elements 
of products of quark fields at finite time- and space separations, their computation requires 
evaluation of functional traces of the single-particle energy and momentum operator, which 



is in general a difficult problem. Using the finite-basis method of fl2fl , we formulate a 
reliable and efficient numerical procedure to compute the Pauli-Villars regularized quark 
and antiquark distributions. The method takes into account the contributions of the discrete 
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bound-state level as well as the Dirac continuum of quarks. 

Finally, using the methods developed in this paper, we compute the isosinglet unpo- 
larized and isovector polarized distribution functions and discuss the results. In principle, 
the resulting distributions should be taken as the starting point for perturbative evolution 
and be compared with structure function data at large q 2 . The evolution of the calculated 
distributions and the comparison with the data will be the subject of a separate investiga- 
tion. Here, we restrict ourselves to a comparison of the calculated distributions with the 
parametrizations of the data at a low normalization point by Gliick, Reya et al. 0, ffl. 



2 Quark distribution functions in the large— N c limit 
2.1 The nucleon in the effective chiral theory 

In the large- N c limit, QCD becomes equivalent to an effective theory of mesons, with baryons 



emerging as solitonic excitations [O, 14]. At low energies, the main guiding principle for 
formulating this effective theory is the dynamical breaking of chiral symmetry, which, in 
particular, results in the appearance of pions as Goldstone bosons. In the long- wavelength 
limit, the effective theory can be expressed in the form of the chiral lagrangian of the pion 
field, whose structure is basically determined by chiral symmetry. The minimal chirally 
invariant interaction of quarks with Goldstone bosons is described by the functional integral 

B0,§ 



exp (iS eS [ir(x)\) = J DipDifj exp (i J d 4 xip(i^ - MU 15 )^ . (2.1) 

Here, ir(x) is the pion field, 

U{x) = exp \m a (x)T a ] , (2.2) 
U^(x) = exp [z7r a (x)r a 75 ] = l±p-U{x) + ^—^U\x). (2.3) 

The quark field possesses a dynamical mass, M, due to chiral symmetry breaking. It is un- 
derstood that, generally, this theory of massive quarks is valid up to an UV cutoff, A 3> M. 
The effective action eq.( |2.3| ) can be derived from the instanton vacuum, where the cutoff is 
determined by the inverse instanton size, and the dynamical quark mass is momentum depen- 
dent. In practice, rather than working with an explicitly momentum-dependent quark mass, 
one usually takes a constant quark mass and applies an UV cutoff to divergent quantities 
derived from eq.( |2.1| ), using some regularization scheme. 

In the effective chiral theory defined by eq. ( |2.1| ), the nucleon is in the large- N c limit 
described by a static classical pion field. In the nucleon rest frame it is of "hedgehog" form 



U c (x) = exp [z(n • t)P(t)] , r = |x|, n = — . (2-4) 

Here, P(r) is the profile function, with P(0) = — tt and P(r) — > for r — > oo, which 
is determined by minimizing the static energy of the pion field. Quarks are described by 
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means of one-particle wave functions, to be found from the Dirac equation in the external 
pion field, 



(*7^-MfP 5 )$„(x,t) = 0, *n(x,*) = expH£J)$ n (x). (2.5) 

which can be written in Hamiltonian form, 

H$ n = E n $ n , H = - J7 V4 + M 7 ^. (2.6) 

The spectrum of the one-particle Hamiltonian, H, contains a discrete bound-state level. 
This level must be occupied by iV c quarks to have a state of unit baryon number. The 
nucleon mass is given by the minimum of the bound-state energy and the aggregate energy 
of the negative Dirac continuum, the energy of the free Dirac continuum subtracted ||, 

M N = N c E lev + N C J2 (En-E^). (2.7) 

E n <0 

(E^ denotes the energy levels of the vacuum Hamiltonian given by eq.( |2.6j ) with U = 1.) 
It is understood that eq. fl2.7Q is made finite by some regularization method, to be discussed 
below. 

Nucleon states of definite spin-isospin and 3-momentum are obtained by quantizing the 
rotational and translational zero modes of the soliton (their contributions to the energy 
are 0(1/N C )), by integrating over the corresponding collective coordinates with appropriate 
wave functions || [F|. 



2.2 Quark distribution functions in the effective chiral theory 

The simplest way to determine the quark distributions inside the nucleon is to use the infinite 



momentum frame and calculate the number of partons there ||11|| . It is well known that the 
infinite momentum helps to separate quarks belonging to the nucleon from the vacuum ones, 
provided that the transverse momenta of the particles are not growing with the nucleon 
momentum ||11|| . In our chiral quark soliton model this condition is, of course, satisfied. 

More precisely, the quark distributions as functions of the Bjorken variable, x, are, by 
definition, the number of (anti-) quarks whose momentum, say, in the z direction, is a 
fraction x of the nucleon momentum, Pjy, in the infinite momentum frame where 

M N v , , 

P N = t=^= , v -> 1. (2.8) 



2 



The number of (anti-) quarks can be expressed through the nucleon matrix element of the 
creation and annihilation operators, a + , a (for quarks), and b + , b (for antiquarks). We define 
the quark and antiquark distribution functions as 

Di{x) = J ^2KS\x- — UNM^)^)\N^ t (2.9) 

Diix) = J^^5(x-^yN v \bi(kMk)\N v ). (2.10) 
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Here % denotes the set of quantum numbers characterizing the quark, such as flavor and 
polarization. 

In the infinite momentum frame it is possible to express these matrix elements in terms 
of the quark field operator, 

^(x,t)= J j^j^J^Yl [a c (k)exp(-zA; ■ x)tt c (k) + 6+(k)exp(i/c ■ ac)v c (k)] , (2.11) 

where u(p), v(p) are the wave function of the free quarks and antiquarks, normalized to 
uu = —vv = 2M, and c = 1,2 denotes the polarization. Let us consider the Fourier 
transform of the equal-time product of ip and tf) + , 

/A 1 ^, ! exp[- ! k.(x I -x ! )]^(x 2 , t )V-(x I , t ) 
= £^[<> c + (kK'(kK(k)Mk) 

c,c' 

+a+(k)6+(-k)<(k)w c /(-k) exp(2z£; t) + 6 c (-k)a c /(k)^(-k)« c (k) exp(-2z£; t) 
+ 6 c (-k)6+(-k)<(-k)^(-k)] . (2.12) 

Averaging this operator over the nucleon state in the infinite momentum frame, with k 3 = 
— v 2 , v — > 1, we get zero for all terms on the r.h.s. but the first one. Indeed, 
the probability to find a correlated quark-antiquark pair with very large opposite momenta 
(the second and third terms) in a fast moving nucleon goes to zero as v — > 1. Similarly, 
the probability to find antiquarks moving in the opposite direction to the nucleon with large 
longitudinal momenta goes to zero as v — > 1 (the fourth term). To be more precise, these 
matrix elements decrease with the nucleon momentum and can, in principle, contribute to 
the structure functions of non-leading twists. However, we are interested in the leading-twist 
distribution functions, and can therefore neglect all terms in eq.( |2.12j ) except the first one. 
The distribution functions can thus be expressed in terms of equal-time products of the field 
operators as 

Dfa) = J -j^yS - JnJ / ^1^2 exp [-ik ■ (x x - x 2 )] 

x(iy v |^(x 2 ,t)r^( Xl ,t)|iv v ), (2.13) 

/($ k ( Jv^ \ f 
(27r)3^ V _ P~ J J d3o ° ld3x2 exp t _ ^ k ' ( Xl ~ X2 ^ 

x(iV v |Tr [r^(x 2 ,t)V(x 1 ,t)] |JV V ). (2.14) 

The flavor and spin matrices, r i; depend on the particular distribution one is interested 
in. For example, for the number of partons polarized along or against the direction of the 
nucleon velocity one should use 

r* = 7°^, (2-15) 

respectively. 

In order to evaluate the equal-time matrix elements in the nucleon state in eqs. fl2.13l , 
2.14|) we consider the more general matrix element of the time-ordered product of quark 
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fields. In the effective chiral theory, this nucleon matrix element can directly be computed 
as a functional integral over the quark and pion field with the effective action, eq. ( |2.1| ). At 
N c — > oo, the integral can be performed using the saddle point method ||. One finds that 



i(iV v |T {V>(xi,*i)^(x 2 ,* 2 )} |iV v ) = G F (x h x 2 ) 



(2.16) 



where Gp(xi,x 2 ) is the Feynman Green function in the background pion field corresponding 
to the moving nucleon. This saddle point solution of the classical equations of motion can 
easily be constructed. Indeed, since the saddle point equations are relativistically invariant, 
it is evident that the pion field is of the form 



vt 



(2.17) 



where U c (x) is the stationary hedgehog pion mean field in the nucleon rest frame. The 
Feynman Green function, eq.( 2.16 ), is thus determined as the solution of the inhomogeneous 
Dirac equation 

d 



dxi 



1 + 75 



f/ 75 (x, t) 



GF(xi,ti,X 2 ,t 2 ) 

r,(x./i • 



= 5^(x 1 -x 2 ), 
f>t(x,t). 



(2.18) 



2 . 2 (2.19) 

With eq. (|2.16|) , the distribution functions, eqs.( |2.13| , 2.14J) , can now be written in terms of 
the Feynman Green function as limits at t 2 — > t% ± 0, 



DAx) 



Di(x) 



-iNr. 



d 3 k 



8 x 



Pn 



Jd 3 x 1 d 3 x 2 exp [— ik ■ (x x - x 2 )] Tr [r^^ (x x , t u x 2 , t 2 )] 



*2=tl+0 



(2.20) 



d 3 k 



5 x 



Pn. 



x Jd 3 xxd 3 x 2 exp [-ik ■ (x x - x 2 )] Tr [TiGp^x, t u x 2 , t 2 )] tl=sh+0 ■ (2.21) 

To compute the quark and antiquark distribution functions one needs an explicit rep- 
resentation of the Feynman Green function in the background pion field corresponding to 
the fast-moving nucleon. We now want to demonstrate that this Green function can be ex- 
pressed in terms of the single-particle quark wave functions, 3> n (x), and energy eigenvalues, 
E n , in the nucleon rest frame, eq. (|2.5| ). The quark eigenfunctions in the time-dependent 
pion field, eq. (12.171) , can be obtained from the ones in the static pion field in the rest frame 
by a Lorentz transformation. We can thus write a representation of the Feynman Green 
function as 

(j?f(xi, ti, x 2 , t 2 ) 

= -iswleih-h) ]T <Mx;)$n(x 2 )ex P H£ n (t;-i 2 )]- 

I noil— occup. 

-0(t 2 - h) J2 <Mx' 1 )<i n (x 2 ) exp [-iE n (t\ - t' 2 )]\ (2.22) 

occup. ) 
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Here t' and x' are the Lorentz transforms of the coordinates, 

, Xi >2 - v£l, 2 £l,2 - V ■ Xi 



X 



1,2 " VT^ ' ^ ~ 'VT^ ' (2,23) 

and 5[v] is the Lorentz transformation matrix acting on the quark spinor indices, 

S[v] = exp Qcro3^) , <V = ^[7m>7i*]> tanh(^) = u. (2.24) 

In the first term of eq.( |2.22| ) the summation goes over non-occupied states, that is over the 
positive-energy Dirac continuum, in the second term over occupied states, that is, over the 
negative continuum and the discrete bound-state level. 

Let us prove that the Green function defined by eq. (|2.22|) indeed satisfies the inhomoge- 
neous Dirac equation, eq. (|2.18f) . We first note that the Lorentz transformed single-particle 
wave functions, 5'[v]$ n (x / ) exp(— iE n t'), satisfy the homogeneous Dirac equation, eq. ( |2.18| ) 
with r.h.s. equal to zero, therefore so does the Green function, eq.( |2.22D , at t\ ^ t 2 . At 



t\ = t 2 the Green function has a discontinuity. Taking the time derivative of eq.( [2.22j ) at 
ti = t 2 we find that the l.h.s. of eq. Q2.18| ) is equal to F(xi, x 2 )5(ii — t 2 ), where 

F(x 1; x 2 ) = 5>xp [-z J B„v(x , 1 - x 2 )] 7°5$ n (x' 1 )<l n (x 2 )5- 1 . (2.25) 

all 

This function should be equal to 5( 3 )(xi — x 2 ) in order to satisfy the inhomogeneous Dirac 
equation, eq.( [2.18|) , for the Green function^. To prove this, we first convince ourselves that 



the sum, eq. (|2.25|) , vanishes for xi ^ x 2 , and then show that it is indeed a delta function 
with coefficient unity. 

We introduce temporarily two moments of "time", t\ = vx' x and t 2 = vx 2 , and note that 
at ti > i 2 eq. (|2.25| ) can be written via the "retarded" Green function, 

F(x 1; x 2 ) = z 7 5G r et(x / 1 ,(vx' 1 ),x 2 ,(vx 2 ))5- 1 , 
G re t(xi, t h x' 2 , t 2 ) = -i9(t! - t 2 ) E $n(x' 2 )$ n (x / 1 ) exp [-iE n {h - t 2 )] . (2.26) 

all 

We need this function in the space-like region, since (ti — t 2 ) 2 — (x' x — x' 2 ) 2 = (v • (x' x — 
x' 2 )) 2 — (x' x — x' 2 ) 2 < 0, at least when xi ^ x 2 . However, the retarded Green function is zero 
in the space-like region. This is obvious from physical considerations: the retarded Green 
function determines the evolution of a wave packet which at t\ = t 2 is localized at x' x = x 2 , 
which cannot reach the space-like region. More formally, this can also be proved using the 
perturbation expansion for the retarded Green function in the external pion field, 

Grcti^-i, h, x' 2 , t 2 ) = G^ ct {x[, ti, x' 2 , i 2 ) 
+ jdV 3 c^ t (x;,^ + ... (2.27) 

The free retarded Green function, G® et , is zero for space-like separations, therefore the second 
term is non-zero only for t± — t% > \x.[ — x' 2 | and for t 3 — t 2 > |x 2 — x 3 |. This implies that 

1 In field theory, eq. (|2.25|) represents the equal-time anticommutator, {ip(xi,t), ^ + (x2,i)} = <5' 3 ' (xi — X2). 
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it is non-zero only for i% — ti > \x[ — x(j| + |x(j — x 2 | > |x^ — x' 2 |. This argument is easily 
generalized to any term of the perturbation expansion. 

Thus, .F(xi, X2) is zero for xi 7^ X2. It therefore must be proportional to a delta function 
in X! — x 2 or its derivatives. This "point" singularity can come only from states with large 
energies, E n . For such states, one can neglect the pion field and replace the wave functions, 
$„(x), by the eigenfunction of the free Hamiltonian. Saturating the sum in the r.h.s. of 
eq. (|2.26| ) by plane waves we obtain 

i ? (xi,x 2 ) 



xi - x 2 ; 



(2.28) 



(In other words, the leading short-distance singularity of the Green function in the back- 
ground pion field is the same as that of the free Green function.) This completes the proof 
that eq.( |2.22| ) is a representation of the Feynman Green function in the time-dependent 
background field, eq. (|2.18|) . 



We now use the representation, eq. ( |2.22| ), to express the distribution functions directly 
through the quark wave functions in the nucleon rest frame. Passing to the Fourier transforms 
of the quark wave functions, integrating in eqs.( |2.2fl| , |2.21| ) over X12, and taking the limit 
v — > 1, we finally obtain 



— — $t (k)(l + 7 7 3 )7°r*^ 3 + E n - xM„)$„(k), (2.29) 

.KU1|.. K^ 71 ) 



/d k" 
— $t (k)(l + 7°7 3 )7 rV(^ 3 + E n + xM N )<S> n (k). 



non— occup. 



(27T) 



(2.30) 



These formulae represent the quark distributions as sums over occupied states, the anti- 
quark distributions as sums over non-occupied states. We can also write an alternative 
representation for the distribution functions, using the time-ordering opposite to the one 



in eqs.( |2.20| , |2.21| ), and the fact that the discontinuity of the Feynman Green function is a 
space 5-function. In this case we get 



DAx) 



—N C M 



E / ?S $ n(k)(l + 7°7 3 )7°r^(fc 3 + E n - xM N 



N 



non— occup. 



(27T) 



(2.31) 



/d k" 
— ^^(k)(l + 7°7 3 )7°r i 5(£ ; 3 + E n + xM N )$ n (k). (2.32) 



occup. 



(27T) 



In both eqs. (|2.29| , |2.30| ) and eqs.( |2.31| , [2.32j ) it is understood that the contribution of free 
quarks is to be subtracted. In the original representation of the distribution functions 
through the Feynman Green function, eqs.(|2.20|, |2.21|), this means that the leading short- 



distance singularity of the Green function in the background pion field is canceled by the one 
of the free Green function, cf. eq.( |2.28| ). However, eqs.( |2.29| , |2.30| ) and eqs.( |2.31| , p.32j ) still 
contain ultraviolet divergences and have to be made finite by some regularization scheme, 
to be discussed in Section 3. 

It is clear that the equivalence of the representations of the distribution functions as sums 
over occupied and non-occupied states, eqs.( p.29| , |2.30| ) and eqs.( |2.3l| , |2.32|) , is based on the 
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completeness of quark states in the external pion field, eq.( |2.2<j| ), which can in principle be 
violated by ultraviolet regularization. Fortunately, it is possible to regularize the theory in 
such a way that this important property is preserved, see Section 3. 

We emphasize that formulae eqs.( [2.29"| , [2.30Q and eqs.( f2.3ij [2.32j ) are identical to those 
which have been derived in ref.|| from a representation of the distribution functions as 
nucleon matrix elements of quark bilinears separated by a light-like distance. The above 
derivation is a demonstration of the equivalence of the two definitions of parton distributions. 

Finally, to get the unpolarized or polarized (anti-) quark distributions corresponding to 
a nucleon state of definite spin and isospin, one has to take the desired combinations of 
the basic expressions eqs. (|2.29 , 2.30|) and average them with the nucleon spin-isospin wave 
function (the procedure is described in 0). In this paper we consider the distributions 
which appear in the leading order of the l/iV c -expansion, the isosinglet unpolarized and the 
isovector polarized one. To obtain the isosinglet unpolarized quark distribution we sum up 
the two polarizations in eq. ( |2.29|) and average over flavor. One finds [[J 



[u(x) + d(x)} 



occup. 



N C M N E 



n 

occup. 



d 3 k 



$t (k)(l + 7° 7 3 )c^ 3 + E n - xM N )$ n {k). 



(2.33) 



From eq.( [2.32 ) one sees that the corresponding antiquark distribution can also be written as 
a sum over occupied states; it is given by the negative of the r.h.s. of eq. (|2.33| ) at x —>■ —x. 
Henceforth, we shall consider eq. (|2.33| ) as a function defined for both positive and negative 
x, and understand that at negative x it describes minus the antiquark distribution. Alter- 
natively, we can use eqs.( [2.30| , |2.31|) to obtain a representation of the isosinglet unpolarized 
distribution as 



[u(x) + d(x)} 



non— occup. 



—p$t (k)(l + 7 °7 3 )^ 3 + E n - xM N )$ n (k), (2.34) 



non— occup. 



which, again, gives minus the antiquark distribution at negative x. In eq. (|2.33|) vacuum 
subtraction is understood for x < 0, in eq.( |2.34j ) for x > 0. 

The analogous expression for the isovector polarized distribution is 



[Au(ie) - Ad(x)} occnp 
1 



occup. 



d 3 k 



$t ( k )r 3 (l + 7 °7 3 )75 ^ 3 + E n + xM N )$ n (k), 



(2.35) 



with the corresponding antiquark distribution given by the same expression with x — > —x 
(vacuum subtraction is again understood). Here, 2T 3 = ±1 for proton and neutron, re- 
spectively. The alternative representation as a sum over non-occupied states analogous to 
eq. p.34j ) can easily be written down. 
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The isovector unpolarized and isosinglet polarized quark distributions vanish in the lead- 
ing order of the l/iV c -expansion. They are non-zero only after considering rotational cor- 
rections, i.e., expanding to first order in the soliton angular velocity, which is O(l/N ), and 
are expressed as double sums over single-particle levels. We shall not consider them in this 
paper. We note, however, that the techniques developed in Sections 3 and 4 can readily be 
generalized to analyze also these "small" distributions. 



3 Ultraviolet divergences and regularization 

The expressions for the quark distribution functions derived in the previous section are ultra- 
violet divergent and require regularization. To be able to compute the distribution functions 
using the effective chiral theory we must ensure that the ultraviolet regularization does not 
lead to violation of any of their fundamental properties. We want to show now that reg- 
ularization by a Pauli-Villars subtraction, which preserves the completeness of the quark 
single-particle states in the chiral soliton, leads to regularized quark and antiquark distri- 
butions satisfying all general requirements. It was noted in Section 2 that the equivalence 
of the representations of the distribution functions as sums over occupied and non-occupied 
states relies on the completeness of the single-particle states. We shall see below that, with 
Pauli-Villars regularization, this equivalence is preserved for the regularized distributions. 
On the other hand, regularizations based on an energy cutoff not only destroy the equiva- 
lence of summation over occupied and non-occupied states, but lead also to other, related, 
unphysical features. 



The divergent contribution in eq.( |2.33j ) comes from the eigenstates of the Dirac Hamilto- 
nian with large energy, \E n \. One may think therefore that a natural way to regularize this 
divergence is simply to cut the contributions of states with \E n \ larger than some cutoff, ojq\ 



£lcv+0 

]£° cup . = N C M N J duSp[5(uj-H)5(u-xM N + P 3 )(l + 1 ° 1 3 ) 

-U>o 

-{H^H Q ). (3.1) 



(Here we have written the sum over states, eq. (|2.33|), as an integral over energy, the integrand 



being a functional trace involving the Hamiltonian and single-particle momentum operator, 
p 3 ||. This form is useful for investigating the ultraviolet asymptotics.) In fact, we shall 
see below that this regularization is unphysical, leading to a number of problems which are 
easily cured by turning to the Pauli-Villars regularization. However, two reasons force us to 
devote some time to the distribution function with the energy cutoff, eq.( |3.1| ). First, in order 
to see explicitly that the Pauli-Villars subtraction cancels all non-physical effects, one has to 
understand precisely what is to be canceled. Second, our numerical method for computing 
the distribution function (see Section 4) involves an energy cutoff in the intermediate stages 
of the calculation (taken to infinity at the end), so it is essential to know the asymptotic 
properties of the distributions with energy cutoff. 

The ultraviolet divergence of eq. ( |3.1| ) can be derived using the technique developed in . 
One replaces the 5-function by the imaginary part (discontinuity) of the quark propagator 
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in the background pion field, 



Writing the quark propagator in the form 

1 oo + H uj + ^°(-i^ k d k + MU^) 



uj-H uj 2 -H 2 to 2 + d 2 k - M 2 -iM 1 k {d k U^)' 



(3.3) 



one can expand in derivatives of the soliton field. One finds that only the first term in the 
derivative expansion is divergent. In this way one easily obtains the logarithmic divergence 
of the distribution function, eq. (|3.1p, in the limit of large energy cutoff, 

[u{x) + d( x )r ° ccup , 

~ N C M N M 2 log ^sign(x)^ J ^0(k 3 - M N \x\) Tr [C/(k)[C/(k)]1 

(io -» oo). (3.4) 

(Here 9(y) is the step function.) If the soliton field, U c (x), is smooth in coordinate space, 
then its Fourier transform, U c (k), decays exponentially for large |k|. This means that, as 
a function of |x|, the integral in eq.fl3.4p also decays exponentially at large |x|. One should 
keep in mind, however, that the asymptotic formula, eq.( |3.4j ), is valid only in the parametric 
range x ~ 1/N C <C u /M N and thus does not allow us to draw any conclusions about the 
behavior of the quark distribution function at larger \x\. 

It is therefore interesting to derive the UV behavior of the distribution function with 
energy cutoff at larger values of \x\. An asymptotic expansion can be performed in the 
domain 1/N C <C \x\ ~ ujq/M^ by computing the UV divergences of the moments of the 
distribution function (see Appendix A). The asymptotic behavior is given by 

W*> + *(*)1SW ~ £ (« + ^) - > - H) «(-«) 



X 



N / x- \m N 

rf 3 xTr \d k U{*)d k W{K) 



(l/iV c < \x\ ~uj /M n ). (3.5) 

It exhibits at large negative x a rather slow l/x 2 -decay, up to the point x = —2uq/Mn, 
where it ends with a 5-function peak. One should keep in mind that this 5-function appears 
only in the asymptotic limit x ~ —ujq/Mn — > — oo. At large but finite values of x ~ —ojq/Mn 
the delta function in eq. (|3.5|) approximates a narrow peak whose width is much smaller than 
cuo/Mjv- This "large negative x" behavior of the distribution function is an artifact of the 
energy cutoff. We shall see below that both the l/x 2 -tail and the delta function cancel in 
Pauli-Villars regularization. 

In the previous section we have shown that the quark distribution functions can be 
represented in two equivalent forms as a sum over either occupied or non-occupied quark 
states. In general, regularization violates this equivalence. Let us regularize the sum over 
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non-occupied states by an energy cutoff similar to eq.(|8.f|). 

UJQ 

[u(x) + d(x)]2 n . occup . = ~N C M N J dtuS P [5(u J -H)5(LU-xM N +p 3 )(l + 1 1 3 : 

Blev+O 

~{H^H ). (3.6) 

One can easily show that this sum over non-occupied states has the same logarithmic di- 
vergence for large loq as the sum over occupied states, eq. ( |3.4| ). This means that the differ- 
ence between the two representations of the quark distribution functions remains finite for 
ujq —>■ oo. The question is whether this finite limit is zero or not. The answer is, surprisingly, 
no. Moreover, this limit can be computed analytically, using a technique similar to the one 
described in Appendix A. One finds 

J^oo {K X ) +dO*0]occup. - H X ) + rf (^)Cn-occup.} 

Thus, regularization by an energy cutoff leads to an anomalous difference between summation 
over occupied and non-occupied states even in the infinite-cutoff limit. 

The deeper reason for the artifacts encountered with the energy cutoff is that this regular- 
ization violates the completeness of the set of single-particle quark states. The equivalence of 
the two representations of the distribution functions as sums over occupied and non-occupied 
states relies on the locality of the equal-time anticommutator of quark fields (or, equivalently, 



of the discontinuity of the Feynman Green function, eq.( 2.28 )). Leaving out the contribu- 
tion of high-energy states one is dealing with an incomplete set of quark eigenstates, which 
results in a modification of the ^-function equal-time anticommutator. In other words, one 
violates causality, i.e., the anticommutativity of the quark fields at space-like separations. 
What is remarkable, though, is that cutoff regularization leads to anomaly-type phenomena 
which persist even in the infinite-cutoff limit, cf. eq. (|3.7|) . Furthermore, eq. (|3.5| ) tells us 
that such regularization always leads to unphysical results for the distribution function, no 
matter which representation of the distribution function one adopts. It should be noted that 
the usual proper-time regularization of the determinant is of this type. 

A regularization which preserves the completeness of states is the Pauli-Villars regular- 
ization, where one subtracts from the divergent sums a multiple of the corresponding sums 
over eigenstates of the Hamiltonian in which the quark mass, M, has been replaced by a 
regulator mass, Mpy. This mass now plays the role of the physical cutoff of the effective 
theory, which was denoted generically by A in Section 2. The coefficient of the subtraction 
is chosen such as to cancel the logarithmic divergence of the distribution function with the 
energy cutoff, eq.(|3.4|). We thus define 



[u(x) + d(x)]™ up { x ) 



UJQ — >oo 



M 2 



lim [u(x) + d(x)r°(x) ~ T7^H X ) + ^)Cu P .(^ 



t\/[2 rr/ 1 "'WJoccup. 
M IV1 PV 



(3- 



Mr 



and similarly for the sum over non-occupied states. One observes that the unphysical phe- 
nomena associated with the energy cutoff — the negative-x behavior, eq.( p.5|) , and the 
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anomalous difference of summing over occupied and non-occupied states, eq.Q3.7p) — are 
proportional to the quark mass squared, M 2 . Thus, the artifacts of the energy cutoff cancel 
under the Pauli-Villars subtraction, eq. Q3.8|) , as it should be. In particular, the Pauli-Villars 
regularized distributions can now equivalently be computed as sums over occupied or non- 
occupied states. 

In this section we have investigated the consequences of ultraviolet regularization on the 
distribution functions using asymptotic expansion techniques. The physical distributions, 
for finite Pauli-Villars cutoff, can only be computed numerically. Below we shall see that 
the numerically computed distribution functions in Pauli-Villars regularization satisfy all 
general requirements. The distributions decrease rapidly for large \x\ and exhibit the correct 
positivity properties, in full accordance with the results of the asymptotic analysis. 



4 Computation of quark distribution functions 

4.1 Spherically symmetric representation for distribution func- 
tions 

We now develop a method for numerical computation of the Pauli-Villars regularized dis- 
tribution functions. Our general strategy will be as follows. We compute the distribution 
functions as sums over quark levels, eqs. Q2.33| , |2.35| ), for a large but finite energy cutoff. Such 



"intermediate" regularization is necessary in order to have expressions which can be com- 
puted using finite basis methods. The physical distribution functions are then obtained by 
subtracting the corresponding sums with the PV regulator mass, Mpy, according to eq. ( |3.8|) , 
and removing the energy cutoff by extrapolation to infinity. In this way, the energy cutoff 
affects only the intermediate steps of the calculation, not the final result. 

For intermediate regularization we now introduce an energy cutoff in Eqs. Q2.33 , 2.35 ) in 
the form 



[u{x)+d(x))* cup , = N c MnY. (n\(l + lV)S(E n -xM N +p 3 )\n)R(E n ), 



occup. 



(4.1) 



[Au(x) - Ad(x)F = -\{2T,)N c Mn 



3 

x (n\5(E n -xM N +p 3 )T 3 (l + 1 1 3 ) l5 \n}R(E n ) 

n 

occup. 



(4.2) 



We have written the matrix elements between single-particle levels in abstract form, with 
p 3 denoting the ^-component of the single-particle momentum operator. Here, R(E n ) is a 
smooth regulator function with a cutoff, E max . For example, one may employ a Gaussian, 



R(E n ) = exp -— M. (4.3) 
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Alternatively, one may use a Strutinsky (error function) regulator of the kind described in 
ref. ||18|| , which leads to more rapidly converging sums over levels. A corresponding regular- 
ization can be introduced also in the sums over non-occupied states. 

Before evaluating the sums over levels, eqs.(fO, |4 . 2| ) , it is convenient to convert them to 
a more symmetric form. In the derivation of the distribution functions in Section 2, using 
the infinite-momentum frame, it was assumed that the nucleon is moving in the z-direction. 
The orientation of the nucleon velocity is, of course, arbitrary, and the distribution functions 
do not depend on it. We can thus write eqs.(|47 



equivalent ly as 



[u(x) + rf(x)]l up . 



[Au(x) - Arf(x)]l up . 



N C M N £ (n\(l + 7 °v • j)8(E n - xM N + v ■ p)|n> R{E n 



occup. 



(4.4) 



--(2T 3 )N C M N 

x H v • T (! + 7°v ■ -y)l 5 S(E n - xM N + p • v)|n> R(E„ 



occup. 



(4.5) 



where v is an arbitrary 3-dimensional unit vector, v 2 = 1. For the isosinglet unpolarized 
distribution this is immediately obvious; in the case of the isovector polarized distribution, 
eq.(|4.5|), we have made use of the "hedgehog" symmetry of the classical meson field, eq. (|2.4j), 



and the Hamiltonian, eq.( [2.6|) , i.e., the invariance under simultaneous rotations in spin and 
isospin space. We can now pass to a spherically symmetric representation by averaging over 
the orientations of v computing the sum over quark levels. Using the identity 



— / dQ v 5(E n - xM 
An J 



N + V • p) 



1 



2|p| 



|p| - \E n -xM N \), 



(4.6) 



and its generalizations, we rewrite eqs.( [4.4| , [4.5|) in the form 



[u(x) + d(x)} 



R 

occup. 



N c Mn E ( n \ Ail) + ^ (2) 7°P • l\n)R{E n 



n 

occup. 



(4.7) 



[Au(x) - Ad(x)]l up . 
= ~(2T 3 )N C M N J2 (^l^ (3) 75P • r + A^ 7 °p ■ 775 p • r + A^ 7 t • j\n)R(E n ) . 



3 ,., 



occup. 



Here, A^ k \(k = 1,...5), are scalar functions of the magnitude of the single-particle momen- 
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turn operator, |p|, as well as of the level energy, E n , and x, 

h 



A^(\p\,E n , 



X 



2p 



|p| - \E n -xM N \) x 



{E n - xM N ) 



iPl 



(E n - xM> 



IPl 



1 3 {E n - xM, 



2|p| 



12 2 

1 1 (E n - xM N 

2 _ 2 |pp 



Pi 



fc = 1 
2 
3 
4 



(4.9) 



These operator functions of |p| are understood in the usual sense, as functions of the eigen- 
values in a basis where the operator is diagonal. 



4.2 Evaluation in a discrete basis 

The distribution functions, eqs.(4"77|, |4.8| ), are sums of diagonal matrix elements of functions 
of single-particle operators between eigenstates of the Dirac Hamiltonian in the background 
pion field, eq. (|2.6| ). To evaluate them numerically we employ a basis of eigenfunctions of the 
free Dirac Hamiltonian, 

H^i = Effo, H = -i 7 A + M7° (4-10) 

The basis is made discrete by placing the soliton in a 3-dimensional spherical box of finite 
radius, imposing the Kahana-Ripka boundary conditions on the surface ||12|| . The eigenval- 
ues and eigenfunctions of the full Hamiltonian, eq.( |2.6| ), are then determined by numerical 
diagonalization in the discrete basis, 

^ ] Hj,jC n j E n c n {, (4.11) 

3 

Hij = I rf 3 x0l(x)^(x), (4.12) 

J box 

$ n (x) = ^c ni ^(x). (4.13) 

i 

Since the operator |p| is a function of the free Hamiltonian, 



|p| = ^m-M\ (4.14) 
it is diagonal in the basis of H eigenstates, eq. (|4.10| ), and one has 



(i|/(|p|)li> = /(IpI0*v, IPl* = \/(Ei° } ) 2 ~ M\ (4.15) 
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for any function, /(|p|). Using this property one can explicitly evaluate the matrix elements 
between levels in eqs.(|4.7|, |4~8| ), and obtains 



[u(x) + d(x)] R 



occup. 



N C M N £ J2 c *m c nj \A^(\p\ i ,E n ,x)5 ij + A^(\p\ i ,E n ,x) ( T °p • 7 )J R(E n ) 



occup. 



I. J 



(4.16) 



[Au(x) - Arf(x)]f ccup . 

= -\(2T 3 )N C M N £ E<^' [A (3) (|p|„^n,^)(75P ■ T)y 

n »,j 

occup. 

+ AW(|p|,, £ n , x) £( 7 °p ■ 7k(75P ■ r) fej + A (5) (|p| 4 , E n , x) ( 7 °r • 7 ), 



i2(£ n ) 



(4.17) 



Now the y4A fe \ (fc = 1, . . .5) are ordinary functions of the eigenvalues of |p| in basis states, 
|p|j. Here, (. . .)«-,• denote the matrix elements of the corresponding operator between basis 
states. We remind that the corresponding antiquark distributions are given, respectively, by 
minus the r.h.s. of eq.( [4.16| ) at negative x, and that of eq.(^T7|) at negative x, see end of 
Section 2.2. 

We note that the Dirac and isospin structures appearing in eqs. (|4.16 , 4.17 ) are essentially 
the same as those in the sums determining the nucleon mass and isovector axial coupling, 
l7|. To simplify the calculation of the matrix elements one may use that 



(3) 
9 A 



7 °p- 7 = F0-M7 , 



(4.18) 



and that the operator appearing in the first and second term on the r.h.s. of eq. ( }4.17[ ) can 
be expressed as an anticommutator, 



7sP • t = - {757°^ ■ 7 , 7°P • 7} 



(4.19) 



The terms in the sums eqs. (|4.16| , |4.17|) are proportional to a step function depending 
on the level energy, E n , the momenta of the basis states, |p|j, and the Bjorken variable, 
x. The expressions can not directly be used for numerical evaluation in a discrete basis, 
since the result would be a discontinuous function of x. There are, however, ways to convert 
eqs.( |4.16| , |4.17|) to a form suitable for evaluation in a discrete basis. One possibility is to 
apply Gaussian smearing in x to the distribution functions. Let us define 



D 



smeared ( 



X) = 



'Jy/ir J-00 



dx' exp 



— ix — x 



r 



A2' 



D(x'), 



(4.20) 



where D(x) stands for the regularized isosinglet unpolarized or isovector polarized distribu- 
tions. Here, 7 1 is a small but finite number. These "smeared" distribution functions can 
now be calculated using eqs. Q4.l6] , [4. 17| ) with replaced by the corresponding "smeared" 
functions in x, 

A^(\p\ u E n ,x) - A^ s ^ eA (\p\ h E n ,x). (4.21) 
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These are now continuous functions of the level momenta and energies, so one may perform 
the sums over levels in the discrete basis, provided one makes sure that the separation be- 
tween the momentum eigenvalues of the basis states is significantly smaller than the smearing 
width, 7. The level spacing is inversely proportional to the size of the Kahana-Ripka box, so 
it becomes necessary to use rather large boxes to attain small values of 7. In the calculations 
described in this paper we use a value of 7 = 0.1, which requires box sizes > 20M -1 . 

In this way one can compute the smeared distribution functions, eq. (|4.20|) . At values 
of x where the exact distributions are smooth the smeared functions provide an excellent 
approximation to the exact ones. An exception is the isosinglet unpolarized distribution near 
x = 0. The exact distribution has a discontinuity at x = 0, which becomes a smooth crossover 
of width ~ I/7 in the smeared distribution. It is possible to recover the discontinuity 
by "deconvoluting" the numerically computed smeared distributions. Dividing the Fourier 
transform in x of the numerically computed smeared distribution by that of the Gaussian, 
eq. (|4.20| ), one can reconstruct the Fourier transform of the exact distribution for values 
of the argument up to ~ I/7. The exact distribution function itself is then obtained by 
inverse Fourier transformation, incorporating the known asymptotic behavior of the Fourier 
transform for large arguments corresponding to a discontinuity at x — 0. 

We thus compute the mode sums for the smeared distributions, eqs. ( |4.16| , |4.17| ), for a 
number of values of the energy cutoff, typically up to E max ~ 10 M, and also the corre- 
sponding sums with the constituent quark mass replaced by the PV regulator mass, Mpy. 
We then perform the PV subtraction, eq. (|3.8|) , and remove the energy cutoff by numerical 
extrapolation to E m ^ — > 00 pointwise in x. One computes a least-squares fit of the PV 
subtracted sums to a constant plus inverse powers of E max , for each x. The stability of the 
extrapolation can be checked by adding more terms to the fit. 

In Section 3 we investigated the asymptotic behavior of the distribution functions with 
an energy cutoff and noted a number of unphysical features, which are removed by the 
Pauli-Villars subtraction. This can also be seen directly in the numerical calculations. The 
numerically computed distribution functions for finite energy cutoff, eqs. ([OB], |4.17| ), exhibit 
a "tail" at large negative x, which is proportional to M 2 , consistent with the asymptotic 
formula, eq.( |3.5|) . (For summation over non-occupied states, the "tail" occurs at positive 
x.) Moreover, the result for the anomalous difference between summation over occupied and 
non-occupied states, eq. (|3.7|) , is confirmed by numerical calculations. Thus, the numerical 
results fully support the conclusions of Section 3. 

Given the equivalence of summing over occupied and non-occupied states in PV regu- 
larization, one may choose any of the two representations for the numerical calculations. 
In practice, it is convenient to compute the quark distributions by summing over occupied 
states and the antiquark distributions by summing over non-occupied states. In this way, 
no vacuum subtraction is required. Furthermore, these sums exhibit asymptotic behavior in 
the energy cutoff earlier than the respective other representations, making the extrapolation 
to infinite cutoff more stable. 

When computing the isovector polarized distribution one must keep in mind that it is 
defined as the limit of zero pion mass of the distribution computed for finite pion mass, i.e., 
for a soliton profile vanishing exponentially at large radii || . (The same limit is understood 
in the definition of the isovector axial coupling, g A .) We thus must carry out the entire 
calculation described above (that is, summing over quark levels, PV subtraction and extrap- 
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olation to E max — > oo) for M T 7^ 0, and take the limit M n — > by numerical extrapolation 
at the very end[|. 

The Gaussian smearing, eq. (|4.20|) , offers a simple possibility to compute the distribution 
function directly as a function of x. We note, however, that the use of the spherically 
symmetric representation, eqs.( 4.1€ , 4.17|) , is not limited to this method. In fact, performing 
other functional transformations of the expressions eqs.( |4.l6| , |4.17[ ) before summing over 
levels one can obtain prescriptions for evaluating the distribution functions in a variety 
of representations. For example, replacing in eqs.( |4.l6| , |4.17| ) the functions by their 
moments, 



A%X\p\i,E n ) 



dxx r 



(m 



1,2, 



(4.22) 



one obtains a formula for numerical evaluation of the moments of the distribution function. 
(Again, one must compute the sums for finite energy cutoff, perform the Pauli-Villars sub- 
traction and extrapolate to infinite cutoff.) We have computed the lowest moments of the 
distributions in this way (m < 10) and verified that they coincide with the moments of the 
numerically computed distribution functions. 



5 Numerical results and discussion 

In the numerical calculations we use the standard value for the constituent quark mass, 
M = 350 MeV, as derived from the instanton vacuum fl?fl . The value of the PV regulator 
mass, Mpv, is determined by reproducing the experimental value of the pion decay constant, 

d A k M 2 . , T M 2 { d A k M 2 PV 



F = AN I AN - 

C J (2tt) 4 {M 2 + k 2 ) 2 C M 2 V J (2tt) 4 (Mp V + k 2 ) 2 

N C M 2 M 2 V 

= l^ l0g ^ (5 ' 1} 

With F n = 93 MeV one obtains M 2 PV /M 2 = 2.52. For the soliton profile, eq.(|J), we use 
the variational form of ref.0, 

P(r) = -2 arctan (^j , (5.2) 

with ro = 1.0 M _1 , which gives a reasonable description of a varitey of hadronic observables 
of the nucleon. For these parameters, the nucleon mass is found to be Mn = 1150 MeV. 
(The nucleon mass is also computed in PV regularization, subtracting from eq. ( |2.7| ) M 2 /Mp v 
times the corresponding expression for the Hamiltonian with Mpy. The contribution of the 
discrete level is also subtracted.) For calculation of the isovector polarized distribution, we 
introduce a finite pion mass in eq.( |5.2D in the form 



= — 2 arctan 



^(l + M^r) exp(-M^r) 



(5.3) 



2 If one computed the isovector distribution directly for a massless soliton profile, P(r) ~ 1/r 2 for r — > 00, 
one would find a singularity at x — (regulated only by the finite box size). The distribution obtained as a 
limit of a massive profile is non-singular at x = 0. 
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This form has the correct Yukawa tail at large r but is not modified compared to eq. ( |5.2j ) at 
r = 0. The limit M n — >■ is taken at the very end of the calculation. 

The result for the isosinglet unpolarized quark- and antiquark distributions is shown in 
Fig.l. For both distributions we show separately the total result (the sum of the discrete 
level and the negative Dirac continuum) and the contribution of the discrete level. One 
sees that the discrete level contributes to the antiquark distribution with a negative sign. 
(The contribution of the discrete level to the r.h.s. of eq. ( |4. 16|) is continuous at x — 0, 
and the antiquark distribution is just given by the negative of eq.( CTB ) at negative x.) An 
approximation in which only the discrete level is taken into account would thus lead to 



negative antiquark distributions |19| . Positivity of the antiquark distribution is naturally 
restored by including the Dirac continuum. This is clear in the light of the discussion of 
Section 3: Restricting oneself to the contribution of the discrete level one is working with 
an incomplete set of states. Only the sum of all levels (discrete plus Dirac continuum) gives 
the correct realization of the distribution function in the effective theory. 

The result for the isovector polarized quark- and antiquark distributions is displayed 
in Fig. 2 (total results and level contributions). One again observes a sizable contribution 
from the Dirac continuum, which reverses the sign of the level contribution to the antiquark 
distribution. Here, however, contrary to the isosinglet unpolarized distribution, no definite 
sign is required a priori. 

The calculated distributions should in principle be used as input for perturbative evo- 
lution, starting with a scale of the order of the cutoff, Mpy ~ 600 MeV. We stress that 
we are computing the twist-2 parton distributions at a low normalization point, not the 
structure function (cross section) at low q 2 , so a meaningful comparison with the data can 
be performed only after evolution to large q 2 . Alternatively, we may compare our calcula- 
tions with the parametrizations of Gliick, Reya et al. 0, |J. Starting from "valence-like" 
(non-singular) quark-, antiquark and gluon distributions at a normalization point well below 
1 GeV, these authors can fit at large q 2 not only all the data in the large-x region, but also 
the recent small-x data down to x ~ 10~ 4 . We emphasize that the quark- and antiquark 
distributions obtained in our approach are precisely of this "valence-like" form. Moreover, 
the normalization points of the LO and NLO distributions of [Q, || are close to our cutoff, 
so one may perform a preliminary comparison without taking into account evolution. 

Fig. 3 shows the isosinglet unpolarized total distribution (quarks plus antiquarks) together 
with the fits of ||. Our distribution is larger than that of || since their fit includes gluons, 
which carry about 30 percent of the nucleon momentum at this scale. For the variational 
soliton profile, eq. Q5.2p, the second moment of the calculated distribution of quarks plus anti- 
quarks is 0.8. (With a self-consistent solution it would be unity, since the energy momentum 
sum rule follows from the equations of motion for the pion field ||.). 

The isosinglet unpolarized valence quark distribution (quarks minus antiquarks) is com- 
pared in Fig.4. Here we have taken in our calculation Mpy —>■ oo, since this distribution 
function is ultraviolet finite and should not be regularized in order to preserve the baryon 
number sum rule. It is interesting to see how this sum rule is realized in the large-A^ ap- 
proach. In fact, the (unregularized) contribution of each level to the baryon number is N c , as 
can be seen, for instance, by integrating eq . (f4.7|) over both negative and positive x. Since the 
baryon number of the discrete level is N c , the contribution of the negative Dirac continuum 
to the baryon number must be zero in order to satisfy the sum rule for the total distribution. 
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Indeed, one observes that the baryon number of the negative continuum vanishes when the 
cutoff is taken to infinity. However, this does not mean that the Dirac continuum does not 
contribute to the valence quark distribution — just the integral of its contribution is zero. 
In Fig. 4 we show both the total result (discrete level plus continuum) and the continuum 
contribution, which integrates to zero. 

The isovector polarized total distribution (quarks plus antiquarks) is shown in Fig. 5. 
The calculated distribution is systematically smaller than the fit of This is related 

(3) 

to the fact that the isovector axial coupling, g A , which determines the normalization of 
this distribution, is underestimated in the leading order of the 1 /^-expansion (with our 
parameters for the chiral soliton we obtain = 0.9). We note that 1/N C corrections to 
this quantity have been computed [20]; the same techniques could also be applied to the 
distribution functions. In Fig. 6 we show the polarized antiquark distribution, which was 
assumed to be zero in the fit of 0. In our calculation it is obtained non-zero, however, 
significantly smaller than the total distribution of quarks plus antiquarks. 

To summarize, we obtain a reasonable description of the isosinglet unpolarized and isovec- 
tor polarized quark and antiquark distributions. In particular, we find a large antiquark 
distribution at the low normalization point, in agreement with the parametrizations of the 
data. 

In the calculations reported here we have chosen the variational soliton profile, eq. ( |5.2| ). 
with a radius which gives a reasonable overall description of a number of hadronic observ- 
ables, for example the iVA-splitting In principle, the classical pion field describing the 
nucleon should be determined as the minimum of the static energy, i.e., as the self-consistent 
solution of the equations of motion of the pion field. The calculation of parton distributions 
with the self-consistent pion field will be the subject of a separate investigation. 

The exact numerical calculations fully support the approximation used in ||, based on 
the "interpolation formula" for the quark propagator in the background pion field. For both 
unpolarized and polarized distributions the differences to the exact results are of the order 
of 10 percent for the total distributions (quarks plus antiquark), somewhat larger for quarks 
and antiquarks separately. 



6 Conclusions 

In this paper we have completed the first part of the program formulated in ||. We have 
computed the leading quark distribution functions in the l/iVc-expansion, namely the isos- 
inglet unpolarized and isovector polarized, in the effective chiral theory. 

Starting from the original definition of parton distributions as numbers of particles frac- 
tion x of the nucleon momentum in the infinite-momentum frame, we have shown that it 
leads in the large-iV c limit to the same results as the QCD definition of distribution functions 
as matrix elements of bilinears of quark fields on the light cone. The fact that the equivalence 
of the two definitions of distribution functions can be established within the mean-field pic- 
ture at large N c shows the scope of this relativistically covariant, field-theoretical description 
of the nucleon. 

We have observed that, generally speaking, the calculation of quark distribution functions 
puts strong demands on the regularization of the effective theory. A crucial requirement is 
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that it should preserve the completeness of the basis of single-particle quark wave functions 
in which one expands the fermion fields in the mean-field approximation. A regularization by 
subtraction, such as the Pauli-Villars regularization, meets this requirement, while methods 
based on an energy (or other) cutoff violate this completeness, and thus, causality In 
particular, the anomaly observed in the difference of summation over occupied and non- 
occupied states in Section 3 shows that one faces here a truly qualitative difference between 
regularization methods, not simply finite-cutoff effects vanishing in the infinite-cutoff limit. 
By explicit calculation, we have shown that Pauli-Villars regularization leads to quark and 
antiquark distributions satisfying all general requirements. 

As to the numerical calculations, we have presented a general scheme for computing the 
distribution functions in various representations. It is remarkable that the Kahana-Ripka 
method, using a basis of eigenstates of the free Hamiltonian, lends itself so naturally to 
the computation of distribution functions after one has converted them to a spherically 
symmetric form. 

The methods developed here, both analytical and numerical, can readily be generalized 
to compute also the "small" quark distributions in the l/A c -expansion, the isovector un- 
polarized and isosinglet polarized distributions. They are given by double sums over quark 
levels in the background pion field (formulas have been presented in ||). Calculations of 
these distributions are in progress. 

We have found reasonable agreement of our results with the fits of Gluck, Reya et al. 
[§, £§]. Indeed, the large- N c approach to parton distributions formulated in || provides 
justification for the picture of "valence-like" distributions at a low normalization point. The 
fact that at large N c the nucleon is characterized by a classical pion field (or, equivalently, 
a polarized Dirac sea of quarks) naturally explains the large antiquark content at the low 
normalization point. The antiquark distributions obtained in our approach are non-singular 
at small x. We note also that the parametric suppression of the gluon relative to the quark 
distributions, which is implied by the effective chiral theory (see Section 1 and ref.@), seems 
not to be in contradiction with the parametrization of the data at low normalization point 
HJ. A 30 percent momentum fraction of gluons at the low normalization point is consistent 
with the suppression of the gluon distribution by M 2 /A 2 . However, in order to make this 
more quantitative one should develop this approach to a level which allows one to compute 
a non-zero gluon distribution. This can be done in the framework of the instanton vacuum, 
using the methods of 0. 
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A Asymptotics for large energy cutoff 



In this appendix we discuss the asymptotic properties of the distribution function with the 
energy cutoff, which was introduced in Section 3 as a device to control the intermediate steps 
of the calculation. In particular, we show how the large-x asymptotic behavior, eq. ( |3.5| ), 
can be derived. This formula assumes the double limit of large x (in the sense of N c x ^> 1) 
and large energy cutoff (ujq ^> M). It is convenient to analyze this limit in terms of the 
moments of the distribution functions. An explicit expression for the moments is obtained 
by integrating eq . (|3~T| ) over x, 
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The delta function of the Hamiltonian can be represented as the discontinuity of the quark 
propagator, cf. eq. (|3.2|). One obtains a representation of the moments as 



M n° = -^j dio'[M n {-iu' + 0) - M n {-iu' - 0)], 
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The integrand of eq. QA.lj) is ultraviolet finite for any fixed u, and an ultraviolet divergence 
appears in eq. ([A.l|) only in the limit ujq — > oo due the large-cu behavior of the integrand 
(typically a power-like growth, see below). Eq. (|A.3|) , on the other hand, contains ultraviolet 
divergences even for fixed u due to large momenta. The difference is that in eq.( A.lj) we 
have the delta function of the Hamiltonian, which constrains the ultraviolet growth of the 
integrand. We therefore have to introduce an additional regularization of the functional 
trace in eq. ([A.3|) at fixed u. The dependence on this additional cutoff cancels, however, after 
taking the discontinuity in eq. (|A.2|) , and thus does not influence the final result. 

Writing the propagator as in eq.( |3.3| ) one can expand eq.( [A.3| ) in derivatives of the pion 
field, 
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For the calculation of the leading large-x asymptotics it is sufficient to restrict oneself to the 
first two terms in eq. ( |A.4| ). The traces are easily evaluated by inserting plane-wave states, 
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fc = 2. 



The integral over q is divergent, but, as said above, the divergent terms do not contribute 
to the discontinuity in ui. Keeping only those terms that lead to the discontinuity we obtain 

if'(w) = AM 2 N c Mlf n (-i) n+2 u n - 3 sign(Reuj) [ d 3 xTr \d k U ~{*)d k U ] (x 
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Inserting this asymptotic behavior in eq.( |A.2j ) we obtain the leading divergences of the 
moments for Uq —>■ oo, 
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Thus, the moments for n > 3 have power divergences with the energy cutoff, while for n = 2 
the divergence is logarithmic. 

The asymptotic behavior of the moments may be expressed in the form of a function of 
x, assuming that 1/N C <^ \x\ ~ ujq/M^. Computing the moments of the function of eq.(|3.5|) 
one may easily check that it corresponds to the large-c^o behavior of the n > 3 moments, 
eq.flAj). 

To conclude, we have shown that, with an energy cutoff, the moments of the distribu- 
tion function generally have power divergences. These manifest themselves not in a power 
divergence of the distribution function at fixed x, but in the occurrence of a "tail" at large 
negative x (large positive x in the case of summation over non-occupied states). The power 
divergences are artifacts of the energy cutoff and cancel under the Pauli-Villars subtraction, 
eq.( |3.8j ). The physical, Pauli-Villars regularized distribution functions are well-localized in 
x. 
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Figure 1: The isosinglet unpolarized quark- and antiquark distributions. Solid line: quark 
distribution, u(x) + d(x), total result (discrete level plus Dirac continuum); dotted line: con- 
tribution of the discrete level (after PV subtraction) to u(x) + d(x). Dashed line: antiquark 
distribution, u(x) + d(x), total result; dot-dashed line: contribution of the discrete level to 
u(x) + d(x). 
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Figure 2: The isovector polarized quark- and antiquark distributions. Solid line: quark 
distribution, Am(x) — Ad(x), total result (discrete level plus Dirac continuum); dotted line: 
contribution of the discrete level (after PV subtraction) to Au(x) — Ad(x). Dashed line: 
antiquark distribution, Au(x) — Ad(x), total result; dot-dashed line: contribution of the 
discrete level to Au(x) — Ad(x). 
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Figure 3: The isosinglet unpolarized distribution of quarks plus antiquarks, ^x[u(x) + d(x) + 
u(x) + d(x)). Solid line: calculated distribution (total result, cf. Fig.|I|). Points: NLO 
parametrization of ref . || . 
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Figure 4: The isosinglet unpolarized valence quark distribution, ^x[u(x) + d(x) — u(x) — 
d(x)). Solid line: calculated distribution (total result); dashed line: contribution of the 
Dirac continuum. Points: NLO parametrization of ref.||. 
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Figure 5: The isovector polarized distribution of quarks plus antiquarks, ^x[Au(x) — Ad(x) + 
Au(x) — Ad(x)]. Solid line: calculated distribution (total result, cf. Fig.|2|). Points: NLO 
parametrization of ref . [|J . 
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Figure 6: The isovector polarized antiquark distribution, \x[/S.u(x) — Ad(x)]. Solid line: 
calculated distribution (total result, cf. Fig.|]). In the fit of ref.f§] this distribution is assumed 
to be zero. 
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